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This constitutes a status report on the research being performed by 
Purdue University's School of Aeronautics and Astronautics for NASA 
Ames/Dryden, under grant number NAG4-1. The topics of research in this 
program include pilot/vehicle analysis techniques, identification of 
pilot dynamics, and control and display synthesis techniques for 
optimizing aircraft handling qualities. The project activities for the 
period of April, 1985 through June, 1986 will be discussed herein. 

The following grant-related publications have appeared in the 
archive journals, and were presented at technical conferences during the 
reporting period 

1. "Closed-Loop, Pilot/Vehicle Analysis of the Approach and Landing 
Task" (Copy attached in Appendix; A). This paper was presented at 
the AIAA Guid. & Cont. Conference in August, 1985. Paper No. 
85-1851. It is also soon to appear in the Journal of Guidance , 
Control , and Dynamics . 

2. "Normalized Predictive Deconvolution. A Time Series Algorithm For 
Modeling Human Operator Dynamics," by Biezad and Schmidt. Journal 
of Guidance, Control, and Dynamics, Vol. 8, No. 6, Nov. -Dec., 

1985. 

3. "Cooperative Synthesis of Control and Display Augmentation," (copy 
attached as Appendix B) . This paper will be presented at the AIAA 
Guidance and Control Conf. in Aug., 1986. 

The current technical activity is directed at extending and 
validating the active display synthesis procedure, and the pilot/vehicle 
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analysis of the NLR rate-command flight conf iguraions in the landing 
task. After evaluation of the NLR data base, attention will turn to 
developing a pilot/vehicle analysis methodology for the lateral- 


directional axis 



Appendix A 



CLOSED-LOOP, PILOT/VEHICLE ANALYSIS 
OF THE APPROACH AND LANDING TASK 


4 




OW6B4AL 1 m 

OF POOR QUAUH 


M.R. Anderson* 

D.K. Schmidt* 
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Purdue University 
West Lafayette, IN 47907 


Abstract 

Optimal-control-thcoretic modeling and frequency- 
domain analysis is the methodology proposed to 
evaluate analytically the handling qualities of higher- 
order manually controlled dynamic systems. 
Fundamental to the methodology is evaluating the 
interplay between pilot workload and closed-loop 
pilot/vehicle performance and stability robustness. The 
model-based metric for pilot workload is the required 
pilot phase compensation. Pilot/vehicle performance 
and loop stability is then evaluated using frequency- 
domain techniques. When these techniques were 
applied to the flight-test data for thirty-two highly- 
augmented fighter configurations, strong correlation was 
obtained between the analytical and experimental 
results. 


Introduction 

One important tool in handling-qualities research is 
the application of closed-loop analysis techniques to 
expose undesirable dynamic characteristics in the 
combined pilot and vehicle system. The pilot has been 
frequently characterized by servo-analytic techniques in 
the frequency domain. One advantage is the results are 
in a form very useful to, and understood by, the flight 
control designer. Much of this work was based on a 
quasi- linear pilot modeling technique developed and 
reported by McRuer, et alJ 1 !. 

A significant contribution obtained from a similar 
technique was furnished by Neal and Smith! 2 * j n 1970. 
In this work, Neal and Smith were able to correlate 
Cooper-Harper 3 ! subjective pilot ratings with 
frequency-domain characteristics of the pilot/vehicle 
system, as modeled, for a precision pitch- attitude- 
tracking task. A single-loop, servo-analytic pilot 
modeling approach was used in this work. 

However, the frequency-domain pilot model 
discussed thusfar has been somewhat limited to single- 
input, single-output systems. Multi-loop models have 
been implemented using sequential loop closure 
techniques, with some limited success. The difficulty in 
the multi-loop case arises in that assumptions are 
required as to the pilot/vehicle system’s loop structure 
and to the proper form of the pilot’s loop 
compensation. There are also several difficulties in 
characterizing task constraints in the frequency domain 
alone (i.e. required bandwidth). 

In the early 1970’s, Kleinman, Baron, and 
Levisonl 4 ! put forth a pilot modeling technique based on 
a time-domain, optimal-control approach. This 
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approach has the potential to be very adaptable for 
more complex, multi-loop situations or tasks. However, 
much of the research based on this pilot modeling 
technique has been focused at attempts to predict 
human operator and vehicle time responses or more 
specifically, statistical performance. Usually, such 
statistical information is not the most useful to the 
flight control designer. 

Recently, Bacon and Schmidt^ presented an 
integrated optimal-control, frequency-domain (OC/FD) 
approach for pilot/vehicle analysis of the precision 
pitch-attitude control task. When applied to the flight 
test results of Neal and Smith, the optimal- control 
approach was shown, not only to agree extremely well 
with the original results presented by Neal and Smith, 
but also to yield additional information on achieveable 
closed-loop bandwidth in the task. This methodolgy 
also provides a quantitative task definition in the time 
domain. The pitch-attitude-tracking task was still 
modeled as a single-input, single-output, closed-loop 
task, as shown in Fig. 1, but the optimal-control 
approach was used to obtain reasonable analytic 
estimates of the pilot’s adaptable dynamics, P(s), in the 
pilot/vehicle systems. 



Figure 1 Neal-Smith/Bacon-Schmidt 

Pitch Attitude Tracking Task 


The research described herein, and in considerable 
detail in Reference [6j, presents an extension of the 
methodology ' of Bacon and Schmidt to the more 
complex approach, flare, and landing task. It is usually 
understood that this task is multi-loop in nature and 
the difficulty has been in the estimation of the 
appropriate pilot loop closures and dynamic 
compensation. In addition, developing a single model- 
based metric correlating with workload has been 
difficult, given multiple pilot loop closures instead of a 
single, combined pilot/vehicle control loop. The 
fundamental approach here is to use the time-domain, 
optimal-control method to estimate the pilot loop 
closures, and then to perform a frequency-domain 
analysis which yields a single dynamically-equivalent 
loop. Finally, the equivalent single loop representation 
is evaluated to extract pertinent vehicle handling 
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quality information. The vehicle configurations in 
consideration here are those initially reported by R.E. 
Smithl 7 ! in the LAHOS (Landing and Approach of 
Higher Order Systems] study. It is again emphasized 
that what is sought in this analysis is an analytical 
methodology that will expose unnaccaptable vehicle 
dynamics in a piloted task and will at the same time 
aid in identifying and understanding desired vehicle 
characteristics. 

The Experimental Data Base 

The LAHOS report summarizes flight tests 
conducted using the USAF/Calspan NT-33A variable 
stability aircraft to study highly augmented fighter 
aircraft in the landing flight phases under VFR 
conditions. Several of the configurations were flown 
with additional control system dynamics described in 
Table 1. These additional dynamics were used to 
represent possible control system dynamics which are 
usually present in highly-augmented aircraft. A 
summary of the thirty-two configurations selected for 
study herein is included in Table 2. 

Table 1 LAHOS Report Control System Types 


! CoDtro! System Types 

First Order Lag 

1 

TjB + 1 

First Order Lead/Lag 

fjS + 1 
foS + 1 

Second Order Lag 
ft = 0.7 

4 + ^+1 

u>i u>, 

Fourth Order Lag, 
Butterwortb Filter 
to Model Digital Lag 
fj — - 0.93, ^ — 0.38, 

u 3 = u A = 16(rad/sec) 

1 

l4 + 2-» + 1||4 + + >1 


The airframe longitudinal dynamics for the 
configurations are representable in the standard form 

* v (t) = A^x v (t) + b v u(t) (1) 

where the state vector is defined by, = (u,w,0,q], 
and the control u(t) is the elevator deflection, 6 e (tj. 
Using the approximations 

T(t) = 0(t) - w(t) (2) 

and 

tW = h(t) - il g/u 0 

the vehicle responses of interest may be written 
y(t) = C v x v (t) + d v u(t) 

where 

y T = |0.q^,n«) . 
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Table 2 LAHOS Configuration Summary 


j Aircraft Characteristics | 

Config. 

Aircraft Dyn. 

] Additional Dyn. j 

No. 

u.Jrad/sec)/^ 

r^erc) 


u> 3 (rad/sec) 

1-C 


0.2 

0.1 

- 

1-1 

1-3 

1.0/0.74 

- 

0.25 

- 

1-6 


- 

- 

16 

2-A 


0.4 

0.1 

. 

2-C 


0.2 

0.1 

- 

2-1 



. 

. 

2-2 



0.1 

- 

2-3 

2-6 

2.3/0.57 


0.25 

16 

r- 

c* 



. 

12 

2-9 



. 

6 

2-11 



- 

16(4th) 

3-C 


0.2 

0 1 

- 

3-1 



. 

- 

3-2 

3-3 

2.2/0.25 


0.1 

0.25 

• 

8-6 



- 

16 

SB 


0.2 

0.1 

- 




. 

- 

1 

2.0/1.06 


0.5 

16 

SSI: 



- 

12 

4-10 



- 

4 

5-1 



. 

- 

5-3 



0.25 

- 

5-4 

5-5 

3.9/0.54 


mm 

- 

5-6 



■H 

16 

5-11 



• 

16(4th) 



16(.5s+l)(.43s + ll 

O- 1 
6-2 

1.9/0.65 

(.2s + l)(l.ls + l)(s* + 2(4)(.7)s + 16) 

(.5s + 1 )(.43s + 1 )(.06s + 1 ) 

(.2s + l)(l.ls+ l)(.ls + 1) 


As previously mentioned, the different 
configurations flown for the LAHOS study include a 
variety of additional control system dynamics. Each of 
these sets of additional dynamics can be represented by, 

x k (t) = A k x k (t) + b k u p (t) 

u(t) = C k x k (t) + d k u p (t) . (3) 

where u p (t) is the stick force applied by the pilot. 

Critical Task Modeling 

Clearly, for aircraft without direct-lift devices, the 
ability to control pitch-attitude is necessary in any 
longitudinal task. But to what extent the pilot has to 
precisely control attitude such that he can effectively 
control flight path, for example, is an interesting 
question to be considered. Here we will attempt to shed 
some light on this question by performing an attitude 
analysis as proposed by Bacon and Schmidt. Since that 
analysis procedure is described in Ref.[5], bow it is 
performed will not be repeated here. However, an 
additional “critical task analysis” will be developed, 
and is one which is intended to give additional 
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information on the vehicle characteristics deemed 
important in the landing task. 

In the approach, flare, and landing task, both 
altitude and vertical velocity are clearly of critical 
interest to the pilot. These parameters are both related 
to tbe vehicle fight path, which is controllable by the 
pilot through elevator position commands (in the 
conventional frontside approach). It is hypothesized, 
therefore, that the ability to precisely control flight- 
pat h-angle is a necessary condition to obtain good pilot 
ratings for that particular vehicle dynamic 
configuration in the approach, flare, and landing task . 

The ability to precisely control flight-path angle, in 
an analytical sense, is equivalent to the ability to to 
minimize the deviation of flight-path angle from a 
desired path. A compatible pilot objective, in such a 
case, may then be stated in the form of a quadratic cost 
function as used in the optimal-control modeling 
approach. An appropriate objective function is then, 

T 

J(Up) = E U im '4r/(qTr* + rUp + gUp)dt} (4) 

T— oo 1 0 * * 


where 


7 e = flight path error 


and 


u p = pilot's control input. 


With the pilot’s objective defined, attention must 
be focused on what observations are available to the 
pilot in the landing task. It is generally accepted that 
pitch-attitude information is very important in any 
longitudinal tracking task. For a YFR task as tested in 
the experiment, vehicle responses available to, or sensed 
by tbe pilot are pitch attitude, pitch rate, and plunge 
acceleration. Also, verticle speed, or sink rate is 
observable. Due to the kinematic relationship among 
these variables, it may reasonably be assumed that the 
pilot could dose control loops based on 0 and B as well 
as 7 and 7 information. Finally, deviation from some 
desired- flight path would be important, and “sensed” 
somehow, at least internal to the pilot himself. 

To complete the critical task definition, some 
appropriate characterization of the commanded flight- 
path to be tracked, 7 c (t), is necessary. Note that with 
regard to the precision tracking performance of the 
closed-loop pilot/vehicle system, it makes little 
difference whether the desired flight path is internally 
generated by the pilot, or is some external exogenous 
command like a flight director, for example. The 
important consideration in meeting our analysis 
objectives is not to actually model a desired glide path, 
but to adequately and consistently represent command 
signal frequency content that is important to the pilot 
in the landing task. 1 

With this in mind, a pre-filter driven by “white” 
noise is used to generate a random pitch-attitude signal 
with frequency content similar to that in the original 
discrete instrument pitch-tracking task used by Neal 
and Smith. This signal is not displayed to the pilot, 
but is used to generate a commanded flight path using 
the vehicle’s 7 /9 response relationship. 


7c(s) _ 1 

*c( s ) T«S + 1 ' 


( 5 ) 
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A value of l/r, 2 = 0.5 sec" 1 was chosen here. In a 
state-space representation, the command signals are 
then cxpressahlc as 


X 


0 

1 0 

X 


0 

\ 

— 

-0.25 

-0.5 0 

K 

+ 

0.25 



0.5 

0 -0.5 

it 


0 


or 

*c( l ) = A^t) + e c w(t) (6) 

where w is a zero mean, Gaussian white noise process. 

Aggregating the command, vehicle, and control 
system dynamics into one state space representation 
gives, 


x v 


A e 0 0 

0 Ay byC^ 

x c 

Xy 

+ 

0 

b ¥ d k 

«P + 

e c 

0 



o 

o 

> 

*k 


b k 


0 


y P ~ I C c C v d v c k! 


x c 

X v 

*kj 


+ [d v d k lu + v. 


( 7 ) 


The command/vehicle system is now in the desired 
form, 


x(t) = Ax(t) + bu p (t) + ew(t) 

y p (t) = Cx(t) + du p (t) + v(t) (8) 

with obvious definitions for A, b, e, C, and d. Finally, 
the pilot’s observation vector, y p is given by 

Yp = be-7e.7,7.M (°) 


where 


7 C = 7 C “ 7- 

These observations lead to the block diagram 
description of the pilot/vehicle system shown in Fig. 2. 
Here the cross-hatched blocks represent pilot 
compensation in his control loops. 

More detailed development and description of the 
methodology can he found in References [5| and |6). It 
will only be necessary to outline a few pertinent points 
here. The optimal control model uses the time-domain 
description of the pilot’s sensing and response 
limitations. Statistical representation of these 
limitations, relevant to this report, are summarized in 
Table 3. The values of r, the pilot’s observation delay, 
and r n , the pilot’s neuromuscular lag time constant, are 
chosen to represent the human operator in his most 
aggressive mode. This is done, coincident with our 
analysis objectives, to expose handling quality “cliffs” 
in the flight configurations by modeling the pilot’s most 
agressive control techniques. 


3 



and 
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C, =(0 d] . 

Also, the A matrix is derived from a Kalman state 
estimator, 


A = A, - E^fV'C, 


where E, is the steady state error covariance and V is 
the pilot's observation noise intensity. 

For the critical task as modeled, the pilot control 
input can be expanded in terms of the pilot’s transfer 
function matrix and pilot observations, 


Figure 2 The Multi-Loop Flight Path Tracking Task 
Table 3 Pilot Model Parameters 


i Pilot Model Parameters 

Parameter 

Value 

Observation Vector 

= I7..T..1.X.M] 

Objective Function 

q„ = (16,0,0,0,0,0) 

Weights 

O 

II 

£ 

*• 

Observation Thresholds 

T-^v = 0*°^ deg 
Tv*# ~ 01 8 

Observation Noise 

-20 dB All 

Ratio 

Observed Variables 

Fractional 

f, = 0.3333 All 

Attention 

Observed Variables 

Observation Delay 

t - 0.2 sec 

Neuromuscular Lag 

r B = 0.1 sec 

Motor Noise Variance 

-25 dB 

Control Input 

Fr (Stick Force in lbs.) 


The coherent part of the pilot’s control response 
can be described in terms of a pilot transfer matrix 
defined by, 


u p (s) = H(s)y p (s) . (10) 

The pilot transfer matrix is obtained directly from the 
optimal control approach solution, 

H(») = 7TT7 K sI " A )/ e (,I_Al,<, d<r(sI-A 1 +b,l e ) 

r n s “ r 1 0 

+sI-A + b,/ e ]' , E,C, T V 1 (11) 


with 


'« = (K 0], 


A, - 


A b 

0 "I/',, 



u p( s ) = [ H i.( s ) +sH - r .( s )b«( s ) 

+ |H 7 (s) + sH 7 (s)Hs) 

+ [H,(s)+sHj(s)]0(s) . (12) 

where, for example, the Laplace Transform definition 

L{m = S0(s) 

has been used for the angular rates. An effective pilot 
transfer function, associated with any pilot observation 
17, can now be defined by, 

P„(s) = H„(s) + sH n (s) 

so that the pilot’s control input becomes a linear 
combination of effective pilot transfer functions and 
observations, 

u p (s) = P^(s) 7e (s) + P 7 (sMs) + P,(s)*(s) . (13) 


The pilot control input above is, of course, 
perfectly compatable with the block diagram 
representation of the flight-path-tracking task 
previously introduced (see Fig. 2). An example of the 
effective pilot frequency responses are shown in Fig. 3 
for LAHOS Configuration 2-2. With the frequency 
responses of the effective pilot transfer functions now 
available, discussion will turn to closed-loop analysis of 
the handling characteristics of the LAHOS 
configurations. 


Analysis and Results 

The actual analysis consists first of an evaluation 
of the attitude dynamics alone, using a variation of the 
Bacon and Schmidtl 5 ! procedure. The intent being that 
obviously undesirable attitude dynamics should be 
exposed at the outset. The hypothesis is that for 
aircraft in which flight-path is controlled by pitch- 
attititude, some definition and evaluation of the 
minimum ability to control attitude is necessary even 
though flight-path control is the ultimate goal. Neal 
and Smith, as well as Bacon and Schmidt, were able to 
analyze the pilot/vehicle handling-quality criteria 
problem in the attitude control task as a trade-off 
between pilot workload required to achieve acceptable 
task performance and a subsequent measure of the 
pilot/vehicle closed-loop performance. Traditionally, 
pilot workload has been shown to correlate with the 
pilot’s phase equalization required. 
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Figure 3 Configuration 2-2 Effective 
Pilot Transfer Functions 

Bacon and Schmidt, as did Neal and Smith, use 
the pilot’s phase compensation, as modeled at their 
closed-loop bandwidth frequency, to quantify the pilot’s 
task workload in the attitude tracking task. They also 
used the magnitude of the resonant peak of the 
closed-loop 0 / 0 c frequency response as a gauge of 
closed-loop performance and stability, in terms of high 
frequency oscillation tendencies of the pilot/vehicle 
system. In both analysis procedures, however, some 
limit was selected on the maximum allowable “low- 
frequency droop” of the closed-loop frequency response. 

More recently, VVaszak and Schmidtl 8 ) have been 
able to characterize the closed-loop performance of the 
pilot/vehicle system through the use of a sensitivity 
parameter defined to evaluate the change in closed-loop 
resonant peak due to a small change (10%) in the 
pilot/vehicle forward path gain. This sensitivity 
parameter, defined by 

SP = (low frequency droop) I 1 (14) 


8 

gave results indicated in Fig. 4 for the original Neal- 
Smith pitch-at titudcvtracking configurations also 
evaluated by Bacon and Schmidt H 


0 








-0 



I© © 'LT& 

1 © . ^ X 


• ao -60 -O -?0 o ?0 40 63 

r 1 LOT PHASE COMPENSATION AT BANDWIDTH (deg) 


Figure 4 -Neal-Smith /Waszak-Schmidt 

Pitch Attitude Tracking Results 

Applying this pitch-attitude-tracking modeling 
methodology to the LAHOS data set reveals the results 
shown in Fig. 5. It is obvious from a comparison of 
Fig. 4 and Fig. 5 that the LAHOS configurations are 
grouped nicely, in terms of Level 1, 2, and 3 ratings, 
but with higher pilot phase compensation required in 
the LAHOS data set (landing task) as compared to the 
Neal-Smith configurations (attitude tracking task). 
This difference is of course consistent with lower short- 
period frequencies corresponding to the lower airspeeds 
encountered in the landing task, thereby requiring the 
pilot to generate more phase lead in landing than in an 
up-and-away flight condition. However, by the 
placement of the pilot-rated Level 1 region in Fig. 5, 
one can conclude that the pilot actually accepts the 
required generation of a modest amount of phase lead 
in the landing task, while the same phase lead 
requirement would result in a Level 2 rating in a pure 
attitude tracking task. It is clear, therefore, that 
although the ratings of the configurations in the landing 
task correlate well with the resonance peak sensitivity 
and phase compensation, the Neal-Smith criteria cannot 
be applied directly to infer or predict ratings in the 
landing task. Adjustments in the allowable phase 
compensation and peak sensitivity (or peak) consistent 
with the results of Fig. 5 is suggested. 



Figure 5 LAHOS Configurations 

Pitch Attitude Tracking Results 
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Attention is now turned toward the analysis of the 
LA1IOS data base in the critical flight-path tracking 
task. In this task, however, directly identifying the 
pilot’s phase compensation is not easily accomplished 
with the current multi-loop block diagram arrangement. 
However, the block diagram in Fig. 2 can be 
manipulated, using the results from Fig. 3, for example, 
to give an equivalent, single-loop representation of the 
flight-pat h-tracking task by defining, 


P.,(s) = P,.(s)G,(s)G,(s) (15) 

where 


G *( s ) “ Q 


and 


= — • 

1“P t( s )G#(s)(-^-(s)) 

f 0 


This equivalent form is now as shown in Fig. 6. 



Figure 6 Flight Path Tracking with 
Equivalent Pilot Function 


Fig. 7 illustrates the equivalent single-loop pilot 
function for LAHOS Configuration 2-2, obtained from 
the vehicle dynamics and from the data in Fig. 3. 
From this figure, one can readily see that the pilot is 
again required to generate lead in the flight-path- 
tracking task, in order to stabilize and maximize the 
tracking performance of the closed-loop, pilot/vehicle 
system. The lead generated (maximum equivalent 
phase peak) is then suggested as a model-based 
indication of pilot workload. 

Bacon and Schmidt also advanced the use of 
required closed-loop bandwidth as an indicator of the 
combined pilot/vehicle system’s ability to track over 
the frequency band. Bandwidth is defined here, as per 
Neal and Smith, as the frequency at which the closed- 
loop phase equals -90 degrees. The closed-loop 
bandwidth, determined from the combined pilot/vehicle 
system, depends, of course, upon both the vehicle and 
pilot dynamics (particularly the neuromotor lag time 
constant, r n ). As the pilot model parameters have been 
held constant for all configurations studied, and the 
smallest achievable r n has been assumed, insufficient 
bandwidth is directly related to sluggish vehicle 
response. Results shown in Fig. 8 verify that low 
comparative bandwidth correlates with higher (worse) 
pilot ratings. However, once sufficient bandwidth is 
attained, other methods must be applied to uncover 
handling quality deficiencies. 



.EQUIVALENT SINGLE PILOT FUNCTION 



Figure 7 Configuration 2-2 Equivalent 
Single Loop Pilot Function 



Figure 8 Pilot Rating and Bandwidth Correlation 
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One of the most important aspects of closed-loop 
performance is the stability robustness of the loop itself. 
Stability robustness, here interpreted as insensitivity to 
small changes in pilot compensation, can best be 
measured using the combined pilot/vehicle open-loop 
Bode, or frequency response. Fig. 0, for example, shows 
the open- and closed-loop frequency responses for 
LAHOS Configuration 2-2. 



-CLOSED LOOP RESPONSE 
.J3PEN LOOP RESPONSE 



Figure 9 Configuration 2-2 Open and Closed 
Loop Frequency Responses 

For good closed-loop stability margins in a 
tracking system, the desired “shape” of the open-loop 
frequency response in the crossover region is well known 

( i.e., constant slope of - 20 dB/decade). Any deviation 
rom the desirable frequency response could then be 
thought of as a reduction in loop quality or 
performance. From Fig. 9, one can see a large peak 
evident in the open-loop frequency response at a 
frequency just greater than the crossover frequency. A 
Nyquist diagram of the same open-loop response (see 
Fig. 101 illustrates that any small phase or gain change 
injected by the pilot could cause an instability in this 
system. The magnitude of this peak as defined in Fig. 
9, herein entitled the high-frequency open-loop peak is 


10 

therefore a measure of loop performance, due to its 
close association with loop stability robustness 
properties. Also shown for comparison in Fig. 10 is the 
Nyquist contour for Configuration 2-1, which received a 
better rating in the task. 



Figure 10 Nyquist Diagram for Configurations 2-1 and 2-2 


The genesis of this peak is significant too. Fig. 11 
shows the frequency response of the LAHOS 
Configuration 2-2 flight-path to stick-force transfer 
function. Considering this figure, and that of the 
equivalent pilot function shown in Fig. 7, one can 
readily see that the high frequency open-loop peak in 
the 7/7 e (s) response (in Fig. 9) is due solely to the 
magnitude distortion arising from the pilot’s frequency 
response. As in any dynamic lead-lag compensator 
design, magnitude disortion of this type will develop 
when the compensator is required to achieve significant 
phase lead. This is easily verified by the large phase- 
lead peak in the equivalent pilot in Fig. 7. Again, this 
phase lead has been introduced in the (pilot’s) 
compensation in an effort to maximize the closed-loop 
bandwidth, and provide stability, within the constraints 
of the pilot’s inherent dynamic limitations. 

A plot of the high-frequency open-loop peak versus 
maximum pilot phase compensation, as in Fig. 12, 
reveals a characteristic grouping of the vehicle 
configurations not unlike those of References (2l,[5) and 
J8]. Those configurations rated best overall (Cooper- 
Harper Level 1) in the approach and landing task are 
appropriately grouped together. The pilot phase 
compensation results for these vehicle configurations 
indicate that a certain amount of effective phase lead is 
acceptable to the pilot in the flight-path-tracking task 
also. Those configurations rated poorly (Level 2,3) in 
the approach and landing task would appear to suffer 
either from insufficient loop quality, as measured by the 
high-frequency open-loop peak, and/or excessive pilot 
workload. 

An additional evaluation of loop quality stems 
from the critical task definition - to minimize flight- 
path error. In fact, the time-domain minimization of 
error has been used in several past efforts as a loop 
performance measure (c.f. Refs. 9-12). 
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Figure 11 Configuration 2-2 Flight Path 

to Stick Force Frequency Response 



Figure 12 High Frequency Open-Loop Peak Results 
for the Flight Path Tracking Task 


11 

This area of performance or loop quality evaluation 
may also be analyzed in the frequency-domain. 
Neglecting noise sources, the steady-state, mean- 
squared flight-path-error may be expressed as, 


a 


2 


= tJ |^M*S»du, (16) 

” o ~e 


where S 7 (u) is the power spectral density of the 
commanded input. One can easily see from Eq. (16) 

above that in order to minimize | — (jo;)|* must 

be small when S 7 (u>) is large. This objective forces 


| (jcj)| 2 to be large (| — (ju>)| 2 « 1 for a tracking 
7e 7 

system) at low frequencies. Therefore, the maximum 
value of the low-frequency pilot/vehicle open-loop 
frequency response will give an indication, in the 
frequency-domain, of the tracking performance. This 
peak value, defined here as the low-frequency open-loop 
peak, can also be plotted along with the previously 
defined workload metric, as in Fig. 13. Here again, the 
characteristic grouping of the thirty-two flight 
configurations is evident, therefore, this low-frequency 
peak is also offered as an indication of loop quality. 
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Figure 13 Low Frequency Open-Loop Peak Results 
for the Flight Path Tracking Task - 


Conclusion 

An optimal-control/frequency-domain (OC/FD) 
methodology is presented, intended as a dynamic 
handling-qualities analysis tool appropriate for such 
complex, multi-loop, man-machine tasks as aircraft 
approach, flare, and landing. When an analysis of just 
the pitch-attitude tracking ability was performed on the 
landing-task data base (L AH OS), a grouping of the 
configurations similar to the Neal-Smith results was 
obtained. However, direct application of the Neal- 
Smith criteria does not appear appropriate for 
estimating the ratings in the landing task since, among 
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other tilings, the pilot appears to accept higher phase 
compensation requirements in the landing task. The 
proposed technique for such an attitude analysis is a 
variation of the method of Bacon and Schmidt that 
does not require any selection of closed-loop “droop" 
nor closed-loop bandwidth to perform the analysis. 

A critical flight-path control task was developed, 
and when a closed- loop analysis of this task was 
performed using the LAHOS data base, excellent 
correlation was obtained with several model-based 
quantities of engineering significance. First, correlation 
was evident between the ratings of the configurations in 
the landing task and the closed-loop bandwidth 
obtained from closed loop analysis of the critical task. 
Also, strong correlation was noted between the ratings 
and the pilot phase compensation, as modeled, and two 
measures of loop quality proposed, all of which result 
from the closed-loop analysis technique. The pilot 
phase compensation associated with a single 
analytically obtained describing function, dynamically 
equivalent to the multi-loop model, was suggested as a 
workload metric. A low-frequency open-loop peak, 
related to the ability to minimize flight path deviations, 
and a high-frequency open-loop peak, indicating 
insufficient stability robustness properties, were, 
recommended as measures of loop performance. 
Finally, the quantities used in the analysis technique 
are fundamental in linear systems analysis, and are 
aimed at providing significant insight into the causes of, 
and remedies for unacceptable man-machine system 
characteristics. 
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ABSTKACT 

The Cooperative Control Synthesis (CCS), 
previously developed to synthesize control 
augmentation so as to optimize pilot opinion 
rating, is extended to the synthesis of display 
augmentation for closed-loop manual control tasks. 
The procedure allows slmulatneous solution for the 
display augmentation and control augmentation 
gains using optimal-control techniques, and 
explicitly includes task related and pilot- 
centered requi rement s in the design objectives. 

Ube of the methodology is demonstrated by 2 

considering a compensatory tracking task and k/a 
controlled element dynamics. Analytical 
evaluation of the various control and display 
augmentation designs synthesized using the 
cooperative methodology indicates that 
simultaneous synthesis may lead to a better 
control display trade-off, ae opposed to designing 
the display after the control augmentation has 
already been designed. 

1. INTRODUCTION 

With the advent of high performance aircraft, 
the amount of information to be processed by the 
pilot to successfully accomplish the assigned task 
has continued to Increase. It has, therefore, 
become more critical to determine the beat 
informational set and display dynamics needed by 
the pilot so as to reduce the pilot ' a workload and 
improve performance. In Ref. {1] it was 
hypothesized that there exists a trade-off between 
levels of control augmentation of a plant and 
display augmentation in terms of pilot workload 
and task performance. That is, a highly 
sophisticated display with minimum stability 
augmentation or a high level of control 
augmentation with lesser display integration or 
intermediate combinations thereof will be equally 
acceptable to the pilot. High degree of 
sophistication in either displays or controls not 
only implies higher cost6 but also the pilot must 
regress further In case of a failure. The 
hypothesis of Ref. [1] was borne out by the 
extensive man-in-the-loop simulation work done by 
12 3] 

Lebacqz et al. * in their effort to quantify 
the di splay/control requirements for helicopter 
approach and landing. They reported that there 
was a di splay/control trade-off for pilot ratings 
of satisfactory, while ratings of adequate-but- 
unsatisfactory depended primarily on the control 
system. Even with flight directors for all three 
axes, the pilots found it difficult to fly without 
adequate flight control augmentation. 


Another area where display control trad*-offa 
may be of Importance Is flight test trajectory 
control (FTTC) • A very high degree of accuracy lm 
following complex trajectories can be achieved ky 
[4] 

automatic controllers , but it is undesirable to 
take the pilot completely out of the control loop. 
Trajectory guidance systems have also been shown 
to lead to discernible Improvements In the ease 
and accuracy with which pilots have approached and 
maintained the desired flight test condition* or 

trajectories^. A combination of trajectory 
guidance and control augmentation may then lead to 
desirable levels of accuracy while still 
maintaining the pilot In the loop. 

Thus there exists a need to synthesize 
stability augmentation and display augmentation 
with special attention to the role of the pilot la 
the system, for a specified task. In thle paper, 
a methodology is presented that Is intended to 
provide a systematic approach to synthesizing 
pilot-optimal control and display augmentation la 
complex, closed-loop manual control tasks. This 
methodology 16 an extension of the cooperative 
control synthesis (CCS) technique developed 

earlier^*^ for control augmentation synthesis. 
The characteristic results from the methodology 
are then Investigated by considering a 
compensatory tracking task beginning with a simple 
2 

k/a plant. Although the methodology la clearly 
Intended for higher order multi-channel tasks, use 
of this plant and task is appropriate to gala 
insight, and for validation studies. 

IK BACKGROUND 

The idea that a control augmentation system 
works in cooperation with the pilot, and a 
technique which considers the augmentation system 
and the pilot to be two controllers working la 

•'parallel' 1 was 6ugge6ted by Schmidt Thia 

cooperative control synthesis technique 
Incorporates a mathematical model of the pilot 
behavior and U6es optimal-control theory to 
synthesize control gains that are pilot-optimal as 
modelled. Since display augmentation, like 
control augmentation, has to be in harmony with 
the pilot's abilities and limitations in order to 
be acceptable to him as an aid in accomplishing 
his task, the cooperative synthesis technique was 
considered to provide an appropriate framework for 
simultaneously synthesizing control and display 
augmentation and thus provide a task specific 
trade-off between control and display 
augmentation. 


* Graduate student, student member A1AA 

**Prof essor. Assoc. Fellow AIAA 


In the Appendix, the problem formulation for 
the extension to the control/dlaplay design caaa 
la presented In its moat general form, and thm 
necessary conditions for optimality are stated. 


1 



Tin* iii'lailb i»l t l*i- diiivutloii can \nt found In Ref • 
|ttl. 'llic itppl i rut 1 on of the methodology to 
|»1 1 nl - 1 1 »- I lit*-- 1 uop synthesis of control and display 
aupjtirnt all on will be discussed in Section 111. 

Other display design techniques based on an 


of the pilot have 

in 


Hess 


and 


I 91 

optimal control model (OCM) * 

been suggested by Lcvlson^^, 

1 12 ) 

Kl elnmnn ct si. All these techniques are 

Intended to lead to flight director designs that 
reduce pilot workload, and are clever applications 
of the OCM. All are Iterative In nature and 
Involve designing the flight director law after 
any control augmentation system has been designed. 
Therefore, these do not explicitly Include the 
display/control trade-off In the synthesis 
procedures themselves. 

An excellent dlscussslon of the functional 
requirements for the design of flight directors 
can be found, for example, In Ref. (13). In 
summary, the fundamental control and guidance 
requirement Is to reduce the tracking error to 
zero In a stable, veil damped and rapid manner. 

The main pilot-centered requirement , whether 
considering control or display augmentation. Is 
that the effective controlled element as perceived 
by the pilot be like k/s over a broad frequency 
region. With such a controlled element, the 
pilot's compensation is a pure gain, which Is 
optimal from the point of minimum pilot 
compensation workload. However, not only must the 
pilot-centered requirements be met, but also good 
overall closed-loop performance must be obtained. 
As will be shown, a performance/workload trade-off 
exists. 

2 Consider a compensatory tracking task with a 
K q / s plant as In Fig.l. As Is well known, such a 

plant Is difficult to control in that It requires 
the pilot to generate lead, which results In 
deterioration in performance and pilot ratings. 
Rate feedback can then be used to provide the 
pilot with lead information, either in the form of 
a "quickened 11 display, or changed plant dynamics 
using control augmentation, or a combination of 
both as shown in Fig. 2. In Fig. 2, K. is the 

G 

display "quickening" gain and is the plant 

augmentation control gain. A feedforward gain 
on the pilot's control Input (6^) is used to 

compensate for the reduced control effectiveness 
of the plant when using feedback control 
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Mg. 1 Compensatory Tracking Task with 1C /•* 
Plant 


With display augmentation only, the effective 
controlled element is : 

K (K.s+l) 

*d <6)/ V e> “ — 

r s 

Thus the controlled element as seen by the pilot 
Is like k/s for frequencies greater than 1/K^. In 
their extensive experimental work, HcRuer et 
nn 

si. 1 J found that for a pure k/s like element, 
the crossover frequency of the piloted system is 
u> 14.3 rade/sec and the pilot's compensa t io^ln 
tfie cross-over region is given by F(s)-K^e 

120.23 seconds for a command input bandwidth of 
1.5 rad6/sec. Thus for this display-augmented 
system to be like k/s at crossover, and be 
etablllzable by the above form of pilot 
compensation, we require that the value of be 
such as to provide a lead of 157 deg. at 
crossover. This leads to the requirement that 

1/K -2.3, i.e. K ~0.4 . The tracking error is 
d o 

then given by : 


with 


e(s)/6 (s) - — =■ 
c / > 


s(s+K K.P(s)) 
o a 


(e +K K J P(s)e+K P(e)) 
o a o 

Thus, for a step command, the error approaches 
zero in the steady-state, which shows that the 
rate-augmented display satisfies the control and 
guidance requirement of Ref. [13]. However, since 
the pilot is not directly perceiving the error 
alone, some loss In tracking performance may 
result. 

Conversely, consider control augmentation 
without an augmented display (K^ c 0). Now with 
rate feedback augmenting the plant, the effective 
controlled element is : 

#w/ V -) “ iOTT T 

r o c 

(here K^l+K^ is assumed in Fig. 2). Thus the 



Fig. 2 Display /Control Augmentation using Rate Feedback 
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t «>4t( i i> 1 ) til v 1 riiiriit in UVi' V/e fur frequencies 
lrt.i. t li/i it K K . Using tin* earn* argument as 
I.cIdU 1 , t In- Vrqul iroii’iit on K for the* system to be 
tiiiililllzaMt by a pure gain ftilot compensation la 


K K Z7.8. 
o c 


Note that these values of K, and K are 

o c 

desirable primarily from the point of view of 
minimizing pilot workload. In reality there would 
be restrict ions on allowable values of these gains 
based on performance requirements and limits on 
loop bandwidth of the control augmented system. 
With display augmentation, if K becomes too large 
the tracking error will increase as the pilot 

would then be making 0 (rather than 0) track 0^ . 

Even though the objective of both the control 
augmentation and display augmentation Is to aid 
the pilot, tlie way tills is achieved Is 
fundamentally different for the two types of 
augmentation as illustrated by the above 
discussion. In display augmentation, the required 
lead is provided to the pilot by a change In the 
numerator dynamics of the effective controlled 
element, while in control augmentation the 
denominator is changed, thus affecting the 
stability properties of the dynamics being 
controlled by the pilot. 


With both control and display augmentation, 
the effective controlled element is : 


x,(s)/6 (b) 
d p 


K K O^s+l) 
o c d 

b(b+K K ) 
o c 


ft may appear as an additional disturbance in fcqn. 

(A.l). 


It is worth mentioning here thot although the 
simplified pilot model Is used in the synthesis 
procedure presented in this report, the complete 
model (with predictor, etc.) is used to evaluate 
candidate designs. Moreover, at each iteration of 
the synthesis process, the parameters (e.g. noise 
intensities) in the simplified mode) are updated 
to yield results that are consistent with the 

complete OCM model. It has been shown^*^ that by 

selecting proper noise intensities Jfor the control 

noise v and the measurement nolee v , 
m y* 

the simplified model may yield the eame human 

operator dynamics as the complete model. In this 

manner the cooperative methodology, as developed 

above, can be used to do "pilot in the loop" 

synthesis of the dleplay/cont rol augmentation 

design. 


IV . AN EXAMPLE 

Consider again the simple compensatory 
tracking t86k, as in Pig* 1, with the controlled 
element (plant) dynamics as considered in Ref.[7]s 

6(s)/fi(s) • K 0 /« 2 (K^-11.7 deg.-sec. 2 /deg.) 

The command 6lgnal 0 Is filtered white noise 
w, with the filter : c 

6 (e)/v(e) - 3.67/(e 2 +3e+2.25) 

c 


Then clearly by letting K Q K C " 1/K^, the 

controlled element can be made like k/s over the 
whole frequency region. But, as stated earlier, 
there would be limits Imposed on values of K and 
K by requirements other than those of just C 
minimizing pilot workload. In this elmple 
example, one i6 already faced with the question of 
how to "best" trade-off control and display 
augmentation. 

111. APPLICATION OF COOPERATIVE SYNTHESIS 


and 

E(w)-0 E[w(t)w(t+o)]-6(o) 


With the tracking error defined as e-0 -8*, It 
1 6 now assumed that the error signal Is displayed 
to the pilot on a "fast" display, with first order 
lag at 50 rads/sec, l.e. 

x d (B) “ 7 b+5o) e<8) 


In the Appendix, a dual performance 
optimization problem Is discussed and the 
necessary conditions for the optimality of the 
various controllers, and the expressions for the 
resulting gain matrices, are derived. The 
interaction between the various controllers is 
shown in Fig. A.l. The association of Controller 
2 (u 2 ) wijth plant augmentation, and of display 
control (u ) with the display augmentation should 
be apparent in this formulation. Moreover, with 

Fj=0 in (A • 1 1 ) and 7^-^^ (^ is the pilot's input) 

alongwith appropriate definition of the system 

matrices A and B. , the structure of the 
o lo 

Controller 1 that evolves from the above 
formulation is similar to that of the Optimal 
Control Model (OCM) - except for the 
simplification that the pilot's observation time 
delay J has been eliminated in the observations 
for Controller 1. The elimination of the time 
delay simplifies the dynamic order of the pilot 
model by eliminating the linear predictor in the 
control structure. The pilot's motor noise (v ) is 
however accounted for In this formulation In fhat 


The pilot's observations are the displayed 

variable x, (in.) (and its associated rate x.) and 
a o 

his objective is to regulate the display signal to 

the best of his abilities. This task is reflected 

as that of minimizing the following performance 

index 



lag time constant, for example, t^- 0.1 seconds* 

For all the analysis resultB presented here, 
the following values were used for the parameters 
in the OCM : 

(1) Pilot is assumed to reconstruct the 

rate information from the dlsplyed variable, 

thus attention allocation is 1 for both the 

observations x, and x.. 

a a 
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(II) otihL-i vtil 1 on thresholds lor * and i 
writ* 11. 01? In. und 0.036 in. /sec. 

i er.prc t 1 Vi- 1 y h.ihcd on visual arc angles of 
0.0!> dig. and 0.1b deg. /nee* at the pilot's 
eye and assuming the display to be 1 ft. away 
f roia t he pilot ' b eye . 

(III) Observation noise ratio of “20 dB for 
each observation. 


performance. The weights Y and Y, 
this ram) control the levels of au 
control energy and "display energy" 
being Introduced into the system. 

J n equal a J . 

2 P 


(scalars for 
ment at Ion 


• respect Ively 


For 


r 2- r 21~°' 


Using the Cooperative methodology, four 
different cases of augmentation were considered. 
These cases are discussed below. 


(Iv) Motor noise ratio of -20 dB. 


A. Pi spl ay Augrocntat 1 on onl j 


(v) Observation time delay, T-0.1 seconds. 

The predicted 'piloted' system performance 
for the unaugmented plant is given In Table I. 
(Note that these results differ from those in Ref. 
|6) because of the different thresholds and 
attention allocation for the display considered). 


Table 1_. Unaugmented System Performance 


8 P 

e rms 
(deg) 

6 rm6 
(deg) 

• 

6 rms J 

(df?g/eec) P 

0.009 

0.8)2 

0.602 

5.93 0.98 


The results in Table 1. are indicative of. 
high levels of tracking error and pilot workload. 
We then wish to improve tracking performance and 
reduce pilot workload by a suitable choice of . 
display and/or control augmentation using only 8 
feedback as shown in the block diagram of Fig. 2. 
The use of only rate feedback is consistent with 
the previous discussion. Note that with display 
augmentation, the signal being observed (and 
regulated) by the pilot is : 

x d (s) “ (s+50) (e " K d® } 
and is no longer Just lagged error. 


Within the framework of the Cooperative 
Display/Control design methodology, the form of 
augmentation shown In Fig. 2 leads to the 
following definitions : 


y i“l x d ,X d^ ^y 


y 2 -lM p ) ; y d -e 

c 2 -[-k c ,k 6 ] ; c d ~K d 


(with as in Eqn. (1)) 


The results of display augmentation 
synthesis, without augmenting plant dynamics, for 
varying F 2rf are listed in Table 11. alongwith the 
evaluation results using the full-order model. As 
la to be expected, the value of K increases as 
F 2<J is decreased, but interestingly reaches m 

limiting value of 0.373 in. /deg. -s ec.,’ 1 as F-.-K). 
Note that this limiting value Is very close to the 
value of obtained earlier for minimum pilot 
workload using the cross-over model. The fact that 
it is slightly lower than the value 0.4 can be 
attributed in part to the lower value of the 
pilot's observation time delay used in the 
synthesis procedure and that error was not a 
primary design objective in the classical 
approach. ThuB as the allowable level of ''display 
energy" into the system is Increased, the 
methodology leads to a design that meets the 
pilot-centered requirements. The results of Table 
II. are also plotted in Fig. 3 and Fig. 4. Fig. 3 
Is a plot of rms error (e) vs. rms manual control 
rate (^), and Fig. 4 is a plot of rms error va. 

total control activity (6) acting on the plant. It 
is clear from these. two figures that the pilot's 
workload (of which 6^ is a measure) and the total 

control activity decrease monotonlcally as the 
display augmentation gain K is Increased. The tub 
error initially decreases aS K is increased, but 
then starts increasing beyond a certain value of 
. The limiting value corresponds to the case 
beyond which increasing will lead to much 
degraded performance without any significant 
reduction in pilot workload. If the pilot's 
performance index is used to predict the pilot 

opinion rating* 11 ’ 17 ^, then the results In Table 
II. also indicate that there will be no further 
improvement in pilot opinion rating for increasing 
beyond the limiting value. 

2L # Control Augmentation Only 


The choice of the cost function should 
reflect the various tradeoffs that are 2 lnvolved in 
the design of the augmentation. For the results 
presented herein, Jf was chosen to be 



with g-.009 (corresponding to for the 

unaugmented case, as in Table I.). Wlerwllle et 

al.^ 16 ^ have shown that the pilot's perception of 
workload is strongly correlated with the rm6 value 
of his control input rate, so the inclusion of 6 
in the cost function J 2 reflects the desire to P 
minimize pilot workload. Clearly the inclusion of 
error reflects the desire to provide good tracking 


The results for control augmentation 
synthesis alone (without display augmentation), 
sre listed in Table III., alongwith the evaluation 
resultB using the full-order 0CM pilot model. 

Again it is noted that as F 2 ~> 0, the limiting 

value of K 10.62 deg. /deg .-sec. that is 

achieved, Corresponds to K K -7.2 which is Quite 

o c 

close to the value obtained earlier using the 
cross-over model with the minimum- pi lot-workload 
objective. Thus as the allowable level of control 
energy into the system is Increased the 
methodology leads to a design that again meets the 
pilot-centered requirements. The results of Table 
III. are also plotted in Figs. 3 and 4. These 
results indicate that the tracking performance 
Improves monotonlcally as the level of control 
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ii iij.im' h K tt I I on lb 1 »,!• iriisnl and ilu* pilot workload 
a 1 t>o ridiKTh. lltiwc Vi* r , tin; total control activity 
art 1 up, on flic plant ( 6 ) bloitu increasing beyond a 
m toll) value of K . A good Intermediate value of 
K oil glit cot i eiiponS to that value beyond which 
veiy high levels of total control deflection 
result, with negligible Improvement In tracking 
performance end no noticeable reduction in pilot 
work! oad . 

Comparing the results for display 
auginentat Ion alone with those for control 
augmentation alone. It Is evident that control 
augmentation has the advantage of significantly 
Improving tracking performance, while display 
augment o 1 1 on has the advantage of significantly 
reducing pilot workload and total control 
deflection. This leads^one to conjecture that 
even for the simple k/s plant, there might he 
some advantage to providing a combination of 
control and display augmentation as opposed to 
control augmentation alone or display augmentation 
alone . 

C. Display Design for Control Augmented Plants 

For four of the control-augmented plants 
above, a display-augmentation synthesis was 
carried out, using the same cost index as 
before. As the level of display augmentation was 
increased (F was decreased) for a given 
control-augmented plant, the same behavior was 
observed as in the case of display design for the 
unaugmented plant, i.e. the error first decreased 
and then started increasing while there was 
continuous reduction in pilot workload (4 ) and 

P 2 

total control activity (6). As with the K Jb 

plant case, the error was found to be the least 
for the choice of F-.* 2 ! for all the control 
augmented plants. The synthesis results and the 
predicted performance corresponding to the "best” 
choice of weighting on the display controller 
(i.e. F «=1) are listed In Table IV. These 
results^are also plotted in Figs. 3 and 4. 

It is interesting to note from Table IV. that 
as the plant itself becomes more and more like k/s 
(or as K is increased), the required level of 
display iugmentation decreases - even though the 
performance index being minimized (J^) Is the same 
in all cases. This is an indication of the 
interdependence between control and display 
augmentation for any given task. Also, the 
synthesized display augmentation Is such that the 
pilot workload (4 ) is roughly the same for all 
V 

the cases presented In Table IV. This result Is 
in agreement with the findings in Ref. [2] that 
the pilot workload tends to be constant for 
specific combinations of control and display 
sophistication. 

D. Simultaneous Control/Display Synthesis 

For the same cost weightings (F^ and F ) 
considered above, the cooperative synthesis^ 
algorithm was again excercised except now the 
simultaneous synthesis of the control and display 
augmentation gains G and G was performed. The 
synthesis results ana the predicted closed-loop 
performance are listed in Table V. Again, as the 
level of control augmentation is Increased, the 


Tab] i* II . |)1 splay Augment f d System 



K d 

e nne 
(deg) 

6 rms 
(deg) 

• 

6 rms 
(dfg/sec) 

J 

P 

10 

0.032 

0.756 

0.749 

5.5B 

0.62 

1 

0.230 

0.623 

0.453 

3.38 

0.41 

0.1 

0.353 

0.634 

0.345 

2.56 

0.35 

0.01 

0.373 

0.638 

0.333 

2.46 

0.34 

0.001 

0.373 

0.638 

0.333 

2.46 

0.34 




Fig. 4 Performance and Total Control Activity 
Comparison 
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X , 

(dt g/deg-sec ) 

K i 

e rmr. 

6 roii. 

4 rmt; 

J 


l 


(deg) 

(dig) 

(dPg/src) 

P 

10 

2 

0.029 

0.116 

0.020 

0.149 

0.700 

0.533 

0.704 

0.622 

5.22 

4.46 

0.70 

0.42 


1 

0.203 

0.283 

0.308 

0.637 

4.4b 

0.35 


0.33 

0.461 

0.841 

0.457 

0.716 

4.36 

0.28 


0.1 

0.617 

1.276 

0.442 

0.774 

4.28 

0.26 


0.01 

0.620 

1.284 

0.442 

0.774 

4.28 

0.26 


Table IV . D1 splay Design for Con t rol - Augmented Plante 

(f m • » 


K 

(deg/de£-sec ) 

K d 

(in/deg-sec *) 

e rmfi 
(deg) 

6 rms 
(deg) 

• 

6 rms 
(dfg/sec) 

J 

P 

0.116 

0.149 

0.123 

0,528 

0.507 

3.61 

0.28 

0.203 

0.285 

0.116 

0.496 

0.525 

3.62 

0.24 

0.461 

0.841 

0,095 

0.454 

0.629 

3.83 

0.21 

0.617 

1.276 

0.087 

0.439 

0.693 

3.86 

0.20 


Table V . Simultaneous Pi splay /Control Design 
<r 2d - » 


F 2 

K 

(deg/deg-eec ) 

K 6 

(in/deg-sec *) 

e rms 

6 rms 

• 

6 rms 

J 

<1 


(deg) 

(deg) 

(dig/sec)-- 

P 

2 

0.082 

0.260 

0.153 

0.547 

0.479 

3.02 

0.29 

1 

0.144 

0.498 

0.132 

0.514 

0.490 

2.76 

0.26 

0.33 

0.344 

1.198 

0.108 

0.469 

0.578 

2.67 

0.22 

0.1 

0.437 

1.576 

0.076 

0.453 

0.647 

2.77 

0.22 


level of display augmentation decreases. The 
results of Table V. are also plotted lo Figs, 3 
and 4. Comparing these results to those of 
Section IV C, above, we note that for the same 
level of closed-loop performance (e rms), the 
total control activity (6) is roughly the same in 
both cases, but the pllot'6 control input rate 

(i ) is lower for the case of simultaneous 
P 

synthesis. This is true even when the results are 
compared for the same level of plant augmentation 
(K ), This can be shown by simply cross-plotting 

thi above results versus K , 

c 

A candidate design that offers improvement in 
all aspects over the unaugmented case Is the one 
obtained by simultaneous display/control synthesis 
with the weights of F »1 and “1 (® ee Table V,). 

From Figs, 3 and 4 ve note that for this case 
(labelled A in the figures), the tracking 
performance is much improved over the unaugmented 


case, the total control activity Ifi much lower, 
the pilot's control input rate is much reduced and 
the level of control augmentation (K ) is not too 
high. The frequency responses of thE effective 

controlled element (x./6 (jo>)) for the candidate 
d p 

design and the unaugmented plant are compared in 
Fig. 5. The candidate design is close to k/e over 
a broad frequency region. (As stated earlier, the 
effective controlled element can be made exactly 
like k/s by increasing the values of the gains E. 
and K, , but doing the former will have the effect 
of increasing the required total control 
deflection, while doing the latter will result in 
degraded tracking performance). The predicted 
pilot describing function (P(Jw)) for the 
candidate design and the unaugmented plant are 
compared in Fig. 6, and the resulting loop 
frequency responses P( Jo>) * x d /^ p ( J w ) are shorn In 

Fig. 7. The cross-over frequency (w^), the gain 


6 


i 




Fifi. 5 Eff ective Controlled Element 

(x d /6 p (j U )) 



Fig. 6 Pilot Describing Fi 6- * Loop Transfer Function 
Function (P(ju)) 



margin (C.M.), and phase margin (P.M.) obtained 
from Fig. 7 are noted In Table VI. We note that 
for the candidate design, not only la the cross- 
over frequency higher, but also the stability 
margins are much improved over the unaugmented 
case. Going back to Fig. 6, we note that for the 
candidate design the pilot's phase compensation is 
much lower than for the unaugmented case. The 
pilot's phase compensation at cross-over, with the 
effect of the neruro-muscular lag and the 

(181 

observation time delay removed 1 J ($ pc) la also 
listed in Table VI. 


Table VI . Frequency-domain Comparison 


Controlled Element 

G.M. 

P.M. 

u 

i pc 


(dB) 

(deg) 

(rads?sec 

) (deg) 

Unaugmented 

6.0 

26 

3.1 

75 

Candidate Design 

9.0 

50 

3.9 

47 


V. CONCLUSIONS 

A methodology was presented that has the 
potential of providing a task-specific, systematic 
trade-off between control and display 
augmentation, taking into account the limits of 
the human pilot as an "inf ormation processor and 
controller". Various control and display 
augmentation laws were synthesized using the 
methodology for a compensatory tracking task with 
a k/s plant. Analytical evaluations of the 
synthesized laws demonstrate the applicability of 
the methodology to meet the pilot-centered 
requirements for display design. The results also 
indicate that simultaneous control/display 
synthesis leads to a better design. as opposed to 


designing the display after the control lava have 
been designed. 

Further theoretical work needs to be done In 
terms of the asymptotic properties of the 
synthesis methodology and the numerical properties 
of the solution algorithm - although no numerical 
problems were encountered. Finally, experimental 
validation is warranted. 
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APPENDIX 

OPTIMAL COOPERATIVE CONTROL/DISPLAY D ESIGN 
METHODOLOGY 

Problem Formulation 

Consider a dynamic system acted upon by two 
controllers, and described by the linear time 
invariant set of first order differential 
equations 

f 

x - A x + B u. + B 0 u 0 + D v (A.l ) 

o 10 A 20 l O 

_ _ m. 1&2 r i 

with xeR n , UjCR , U 2 G R and w a zero-mean 
Gaussian white noise process with intensity W. The 
two controls represent two physically independent 
controllers. 

The display dynamics are assumed to be of the 

form 
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(A. 2) 


wJ i h x rf f K # 
"quickening 
the optimal 
(u and u ) 

law u . 
a 


m Vd " 


B u 
do d 


— d — 

u .f H , where u. Is tl»e display 
" controller. The objective la to find 
cooperative controllers 1 and 2 
along with the optimal display control 


Controller 1 (u ) la assumed to have noisy 
observations available for feedback given by 


y. * C, x 4 C x + C u, + v (A. 3) 

1 lo did ud y 

where v Is also a zero-mean Gaussian white noise 
process y with intensity V . This controller will 
have the form of an LQC Compensator, i.e. It 
consists of full state feedback implemented using 
a Kalman estimator. 


The augmentation controller and the 
display control law u^ ajre assumed to have nolee- 
free measurments y^ and y . , respectively, 
available for feedback, where 



y 2 


C 2o XJ 


F d 



(A.«) 


Fig. A.l Block Diagram for Cooperative 
Control /Display Augmentation 


Note that the above formula tion__does not allow 
feedback of the display £tates x rf to the 
augmentation controller u 2 » Finally, these two 
latter controllers are constrained to have the 
direct output feedback form 

“2 “ V2 “ G 2 C 2o* 

"d “ Vd " G d C d 

The interaction between the different controllers 
is shown in the block diagram of Fig, A,1 

Design Objectives 

Controller 1 is to be optimal with respect to 
the cost 


r— 

x 


(A. 5) 


T 

J l “ T (A,6) 

“l R l V^2 F l“ 2 7dt ^ 

in the presence of the action of control inputs 
and u, . Here E{*} indicates the expected value 
operator and the weighting matrices are > 0, 

Q Jd > 0, Rj > 0, Fj > 0. Conversely, Controller 2 

(u ) and the display control law are to be 
optimal with respect to the cost 

T 

J 2 -E{lim Y / <*\ 0 ^jQ 2d V“l R 2 “l + (A * 7) 

o 


u 2 F 2 u 2 +u d F 2d u d )dt) 

in the presence of the control action Uj« The 
weighting matrices are Q„ > 0, Q« . > 0, > 0, 

F 2 > o. F 2a >c. 2 ° 2d ^ 


Augmenting the system dynamics (A.l) with the 
display dynamics (A. 2), the state-space 
description of this augmented system Is 



Defining x m COL (x, x^), (A. 8) can be written la 
a compact form with appropriate definitions for 
the matrices as 


X - A X + BjUj + B 2 u 2 + B d u d + Dw (A.9) 
The measurmente .can similarly be written as 

y, - c/x + c u . + v 

1 1 u d y 

y 2 m c 2 x ( a . 10) 

F d ■ V 


The two cost functions can then be expressed 
in terms of the augmented state vector x as 

T 

Jj - E{lin i / (X T QjX + u J R j u i + 

T 4 ® o 

T 

J 2 - E{lin I /(x T Q 2 X+^]'r 2 u 1 + (A. 11) 

T* 4 ® o 

“2 F 2“2^ F 2d“d >dt) 

where the weighting matrices Q. and Q are 
appropriately defined. Note that this formulation 
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lb loiiu.il ly i twit lot a will 1 1-player non-iero iub 

119] 

gjtiiu*, and vv seek a Natili solution . 

Sol u t 1 on fo r Uj 

ln_th<* presence of tin* action of control 
Inputs u 2 and u , , bk given by (A. 5), the dynamic* 
of the augmented system (A. 9) are 


where 


X • A x 4 B.u. 4 Dw 
aug 1 1 

y. • C x 4 v 
J 1 aug y 


*.«« - <* * »2 C ! C 2 - Wd 1 

c . ug i (C 1 + Wd' 


(A. 12) 


(A. 13) 


and the performance Index Jj becomes 


T 

Jj - EUim | /(x^Qj-fC^GjFjGjCj)^ + (A.U) 

T-*^° o 

U l R l U l> dt ^ 


T r 


Equations (A ,12) and (A. 14), In the case__of 
uncorrelated process and measurement noises (w and 
v ) and for V >0 (l.e. V - positive definite), 
describe the standard non-singular linear 
quadratic Gaussian regulator problem for 
controller Uj . The optimal ^controller Is 

knovn^ 2 ^ to have the form 


J 1 




(A. 15) 


K i - - r I 1 b I p 


with P >0 and symmetric, the solution of the 
algebraic Rlcatti equation 

*Lg P P PA .„ 8 * «1 +C I C I F 1 C 2 C 2> " <*•”< 

PB^R^BjP - 0 


The dynamics of the Kalman state estimator 


are 


X “A x + B.u. + M,(y t - C x) (A. 18) 
aug 11 1 aug* 

where the Kalman filter gain matrix is given by 

(A.19) 


H 


1 


I C T V -1 
aug y 


with t (> 0) the solution of the algebraic Rlcatti 
equation 


A I + Z A* + DWD T - 
aug aug 

rc T v _1 c r - o 

aug y aug 


(A. 20) 


has the form 


K 1 X; 


Aj X 4 M J y ] 


22 

(A. 21) 


where A — (A +B K -H.C ). Then in the 
1 » aug 1 1 1 aug __ 

presence of this control action Uj , the system 
dynamics (A. 9, A. IS, A. 18) can be written In terms 

of the augmented state vector *q ~ COL (x* X) ss 

A 

(A. 22) 


M i c i 


H.C 
1 u 



which can further be written in a compact form, 
with appropriate definitions of matrices, as 


q - A^q 4 + B^u^ + D'w' 


(A. 23) 

rw oi 

The intensity of the process w' is W' » 1^ v 1* 

The_ index of performance to be minimised by 
u 2 and then becomes 


J 2 - E{lim i /(qVq 4 u^F^ 4- (A.24) 

T+® o 


u d F 2d u d )dt) 


with 


where x 1b the minimum mean-square estimate of the 
system state vector x* The gain matrix K^ is 
given by 

— 1 T 

(A. 16) 


Q' r 


*2 

0 


K I P 2 K 1 


The design objective can then bje stated as that of 
finding the optimal controller u 2 and optimal 
display control which minimize the cost as 

given by (A. 24), subject to (A. 23). 

[81 

It can be shown 1 J that the gains G 2 and 
which correspond to the simultaneous optimality of 

the two controllers u 0 and u, are given by 
* r r t 


rWs) 


and 


‘4 

0]HL 

4 

0 

Uc 2 o] i 

0 

' B d‘ 

T 

4 


4 

H.C 

HL 

0 

(ic d 0) L 

0 

1 u 


_ 




G » -F 
d r 2d 


, T 

Here, L - E{q q ) satisfies the relation 

aT 


A L + LA* + D'W'D' 1 - 0 
c c 


r ! (Aw26) 


(A. 27) 


and H satisfies 


S olution for u 2 and 


The optimal controller Uj as derived above 


A H + HA 4- Q - 0 (A. 28) 

c c 

where the following definitions have been used 


ORIGINAL PAGE fs 
OF POOR QUALITY 



A - 

c - 


Q z Q' ■» 


M.C 

1 aug 


Vl 


C K r 2 C 2 C 2- C K F 2d C d C d ° 


M l*i* Vi* 1 opniriit of a FI 1 k 1 1 1 l** 1 Maneuver 
Autopilot for a Highly Maneuverable 23 

Aiinaft, 11 A1AA Paper 83-0061, Jan. 1983. 

|5] Duke, E.L. , Swann, M.K., fcnevol dson , E.K. , 

and Wolf, T.D., "Experl enre with Flight Teat 
Trajectory Guidance," Journal of Culdance 
and Control, Vol . 6, No. 5, Sept. -Oct. 1983, 
pp. 393-398. 


Tin* solutions (A. 2b) and (A. 26) are derived from 
the gradient conditions 


tt 2 ■ 2 <VV c 2 01l 


I 0) HL 


) - 0 


(A. 29) 


and 


S- - 2 < F 2d G d lC d 01L 

d 


M.C 
1 u 


HL 


) - 0 


(A. 30) 


respectively . 


|6) Schmidt, D.K., "Optimal Flight Control 

Synthesis via Pilot Modeling," Journal of 
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Thus the solution to cooperative 
control/display synthesis problem requires 
simultaneously solving two Algebraic Rlcattl 
equations (Eqn. (A. 17) for the control gains for 
Controller 1 and (A. 20) for the estimator gains), 
two Lyapunov equations (Eqn. (A. 27) for the system 
covariance matrix and (A.26) for the matrix of 
Lagrange multipliers) and two gradient conditions 
(Eqns. (A. 29) and (A. 30)) which are necessary for 
the optimality of the gains G„ and . A computer 
program using a conjugate gradient search 

( 21 ] 

algorithm with cubic interpolation was 
developed to solve for the optimal augmentation 
gains. 
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